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We show that there is an effect of nonextensivity acting upon the BCS model for 
superconductors in the ground state that motivates its study in the Tsalhs' statis- 
tical framework. We show that the weak-coupling limit superconductors are well 
described by g ~ 1, where g is a real parameter which characterizes the degree of 
nonextensivity of the Tsallis' entropy. Nevertheless, small deviations with respect 
to g = 1 provide better agreement when compared with experimental results. To 
illustrate this point, making use of an approximated Fermi function, we show that 
measurements of the specific heat, ultrasonic attenuation and tunneling experiments 
for tin (Sn) are better described with q = 0.99. 
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1 Introduction 

About 10 years ago, Tsallis[l] proposed the following entropic form: 



where fc^ is a positive constant, p is the matrix density of the system and g is a 
real parameter. The thermodynamics constructed from (1) has been applied to 
several physical systems in the last years. Indeed, nonextensive statistics was 
used to investigate physical systems which presented nonextensive features. 
Amongst them, there are stellar polytrops[3], solar neutrino problem[4], galaxy 
clusters[5], etc. It should be remarked that the formalism established above 
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seems to be appropriate particularly to systems interacting via long-range 
interactions. 

Furthermore, formal results were achieved and techniques were developed. 
Indeed, from the properties of the Tsallis' entropy, Plastino and Tsallis[6] 
generalized the Boguljubov inequalities and justified the variational method; 
also, this nonextensive thermodynamics preserves the Legendre Transform for- 
malism[7,ll]. It's also possible to establish a relation between the partition 
function with g = 1, as usually calculated by the standard statistics, and the 
generalized qi^\ nonextensive partition function[8] (for this purpose one use 
the Hilhorst integral ). This result was crucial to the generalization of the 
Green's functions method[9,10], which is a useful tool in the discussion of a 
many-body system for any possible real value for where q measures the 
nonextensivity of the system. 

In this article study the BCS (Bardeen, Cooper and Schrieffer) model in Tsal- 
lis' statistical framework. Indeed, the BCS theory for superconductivity de- 
scribes very well, at least qualitatively, a vast amount of experimental results 
for a wide range of materials; these materials, henceforth, will l)e called usual 
superconductors, since they are well explained by the BCS theory. Assuming 
that the exchange of phonons between electrons is responsible for the super- 
conducting phenomena, we shall also assume that the criterium for the appear- 
ance of superconductivity is that the Coulomb repulsion is surmounted by the 
effective attraction caused by the phononic interaction. Instead of calculating 
the interactions involved by "first principles" (divergences in the perturbative 
chain postponed for years the development of the theory) BCS simply consid- 
ered a mean-field approach with excellent results. Such improved theory will 
not be considered in this article. In fact, our aim is to develop a ^-dependent 
formalism which may be applied to materials or a class of materials which 
systematically deviates from the usual theory. As a first approach, we will 
concern ourselves to the weak-coupled limit of the BCS theory, but the strong 
coupling limit and the generalization of the theory to other potentials with the 
correlations included, for instance, is already under study and will be subject 
of further publication. 

In next section we show that there is an effect of nonextensivity acting upon 
the BCS model at T = 0, then we generalize the well known gap equation 
and use it to calculate quantities that will be compared with the experimental 
results for the usual superconductors. 
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2 Motivation 

Consider the well known "reduced" Hamiltonian proposed by BCS and used 
to study the appearance of superconductivity in simple metals [14] , 



where the operator b]. = clcL^ creates a Cooper pair in the singlet state and 
the model describes interacting electrons forming singlet Cooper pairs. This 
assumption may be justified in terms of an instability of the Fermi sea be- 
ing caused by binding of zero-momentum pairs. The only relevant electrons to 
the superconducting state are those in a small neighborhood of the Fermi level 
within a layer twice the Debye frequency thick (If one consider electron-phonon 
interaction as the sole responsible for the superconductivity phenomena.). No- 
tice that we have used the convention that an explicitly negative k has down 
spin and vice versa. 

The transformed Hamiltonian, in a semi-classical view, has its unperturbed 
Fermi sea characterized by "spins" aligned only in the z direction; once the 
phononic attraction between electrons is present, a non-null component of Sk 
perpendicular to z shows up and this configuration is determined in terms of 
the orientation Sxk/Szk = tanOk- Indeed, following Bogoliubov's notation[16], 
sin 9k = Mfcf fc and this quantity is related to the superconducting gap energy, 
Afe, this gap may be seen as the order parameter of the phase transition 
superconductor-normal {Ak ^ represents superconducting solutions.). 

Following the prescription briefiy introduced above, the system may be ana- 
lyzed not at an arbitrary temperature T below the transition temperature T^, 
but particularly at T = 0, where all relevant electrons are in the condensate 
state, henceforth, paired. When the temperatures increases, the possible exci- 
tations of the system are the quasi-particles and the excited pairs. As pointed 
out by Schrieffer[17], at the zero-temperature the operator rik -l- may be 
replaced by 2bl.bk, twice the occupation number of pairs. Hence, we have to 
work in the grand-canonical ensemble adding to (2) the operator —/iN, where 
N is ^k + TT'-k and the chemical potential, may be regarded as essen- 
tially the Fermi energy [18]. Therefore, we get a modeled system of interacting 
electrons with net energy — /x, what is consistent with the picture of only 
relevant electrons for superconductivity possessing energies in the vicinity of 
the Fermi level. 

In order to study the attractive potential, we perform a Fourier transform 
of (2) assuming an isotropic superconductor. The 1/2 "spin" operators are 
defined by 
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bi — Six + iSiy — ; (3) 
bj — Sjx — iSjy — Sj . (4) 

After writing the "reduced" Hamiltonian at T = in terms of these pseudo 
spins(3, 4), it is possible to show that the model is isomorphous to the XY 
model: 

iired = + Vij)bibj , (5) 

ij 



where 



^^.• = ^E(^^-/^)e"*^""'"^^ (6) 



and 



(7) 



kk' 



with Q, being a normalizing constant. 

As pointed out by Tsallis and Anteneodo [19], this classical Hamiltonian may 
be studied with Vij oc r~j°' as the coupling between sites, where r.^ = Ti—rj and 
a is the range of the interaction essentially. A typical quantity that should be 
calculated at T = is dr r'^^^r^"', which is related to the internal energy per 
particle. If < a < d. It's a trivial task to show that the internal energy per 
particle diverges in the thermodynamic limit, what was called weak violation 
of Boltzmann-Gibbs (BG) statistics (in contrast to "strong" T ^ violation). 
This result also reveals the role played by the range and the dimension of the 
system in the scenario of the nonextensivity. 

Following the path established above and assuming d = 3, we will show that 
there is a nonextensive effect acting upon the BCS model at T = 0. For this 
purpose we will use the well known BCS mean-field interaction 



kk' 



-V , we \ek- \ek' - < hwD 
, otherwise 



where hwj:) = ksOn and Od is the Dcbye temperature. The spatial dependence 
of this interaction will be found approximating the sums over k and k' in (7) 
by volume integrals. The region of integrability, of course, will be constrained 
to values of k related to the relevant electrons only: ao = > | |k| — /c^l. 
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where me is the electron mass, what defines a neighborhood around kp- Hence, 
we finally get 

V^J^^^v{r^), (9) 



with 



1 cos(/,>-/-j) siii(ni)/-j) 

"('■"^^^fl 

ao cos{aD'r/3) coslkpr/s) + kp sm{kFrf}) sin(Q;i3r^) ] (10) 
and P = i,j indistinctly (Notice = |r^|). 

There's a one-to-one correspondence between the i's lattice sites and the pairs 
at T = 0, so it's reasonable to replace sums over i and j by integrals. Again, we 
get a product of identical integrals and the potential energy is proportional to 
[/ d^rvp]"^. In agreement of what is expected for a extensive system, the first 
part of the above integral converges to a constant factor in the thermodynamic 
limit ( since kp-iCiD > and limAr^oo / d^r {cos{kFr) sm{aD'r)/r^) — 7r/2 ). 
However, the second and third contribution give raise to oscillatory fiuctua- 
tions in the potential energy according to the number of pairs in the system: 



J d^r — [ kp sn\[kFr) sm.{aDT) — an cos{aDT) cos{kFr) ] = 
47r { sm[{kF + aD)N] - sm[{kF - aD)N] } . (11) 

These fluctuations arc bounded but still present even for a very large number 
of pairs and they don't decrease as you add electrons to the system. Further, 
the energy cannot be expressed as the energy per particle times the number 
of particles itself, as commonly done for extensive systems. 

Therefore, we conclude that there is a formal nonextensive effect acting upon 
the BCS model at T = as, mentioned above. The nonextensive character 
may be regarded as a motivation for the analysis of the BCS theory in Tsallis' 
statistical framework. 



3 The Method 



In this section we will generalize the BCS gap equation in the Tsallis' statistical 
framework making use of the well known Gorkov's anomalous functions. Also, 
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it's that the lifetime of the quasi-particles is infinity, since the interaction is 
time-independent. 



Originally, BCS proposed a wave function for the ground state, lipo), describ- 
ing N paired electrons. This wave function was given in terms of a variational 
parameter that should be found minimizing the Helmholtz free energy. Re- 
garding BCS, the superconducting gap energy. A, may be considered the order 
parameter of the phase transition. Notice that even before the full develop- 
ment of the theory this gap energy was observed in several experiments such 
as specific heat, ultrasonic attenuation rate, tunneling measurements, etc. The 
self-consistent equation for the superconducting gap was given by BCS [14] in 
a Hartree-Fock approximation and may be written as Aj^ = ~J2kk' ^kk'ibl') , 
where (6|.) is the expectation value in the BG statistics for the temperature 
dependent creation operator of a Cooper pair. A generalized version of this 
quantity, which includes Tsallis statistics and the g-dependent expectation 
value, shall simply be 

A, = -Y,Vkk'{bl,),. (12) 

kk' 



Through a canonical transformation the system may be simply described by 
quasi-particles with a given dispersion relation[16] In order to calculate the 
g-expectation value found in (12) we define , for q = 1 and except for a mul- 
tiplicative factor, the propagator ^^(k, r) = {T[cka{T) cl^iO)]), where T is 
the temporal ordering operator and h = 1 from now on. The propagator can 
be rewritten as Ga(k,T) = 6'(r)( Cfeo.(r) 4^(0)) - 6'(-t)( 4^(0)cfecr(r) ) in a 
fermionic system. In order to calculate G in interacting systems we may use 
two methods mainly: Feynman diagrams and the construction of the equa- 
tion of moment. Wc attempt to generalize the latter by closely following the 
development discussed in Ref[20]. 

Consider the advanced(a) and retarded(r) double-time Green's functions [21,9]: 



G^l'§^ = 9it-t'){[A{t),B{t')]), (13) 



G'2'^> = -e{t-t'){[A{t),B{t%), (14) 

where 9{r) is the characteristic function times an arbitrary multiplying factor. 
By definition, ([A, B]^)^ = ( AB + r]BA)g, with rj = -1; if A and B are 
bosonic operators, and rj = 1 if fermionic. The Green's functions admits a 
Fourier series expansion with Fourier transform 

1 7 

^A,B:^^'^^=- j dTG^j:i{t,f)e^-\ (15) 
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where j = r, a and r = t — t' . Assuming that [A(t), B(t')]^ is a non null 
operator, the Green's function can be represented by means of the Hilhorst 
integration [9, 8], which relates the g-expectation value of observables to its 
standard q — 1 analogous average in the BG statistics. The correlation func- 
tion is given then by: 



B{t')A{t)) 



Z, 



1 



oo 



where the kernel Kq{^) is a ^'-dependent function presented in the integral 
representation of the Green's function, according Lenzi's notation[10]. 

We shall calculate the equation of motion of each propagator of interest regard- 
ing the phenomena of superconductivity and than calculate the gap energy, 
in order to compare it with experimental measurements. 



In a time-independent Hamiltonian (such as the case) the corresponding Green's 

--fcT 



functions depend only upon t—t' = r. The equation of motion for -C Cfe|, cL 3>^-*, 



is: 



where jk = — A*- For <^ c^kiclp the procedure is identical and 

(w + 7fe) < cL,^, 4^ >«= A < 4p c_fct . (18) 



This is the Gorkov's anomalous function and it's strictly related to the su- 
perconducting phenomena, i. e., the appearance of Cooper pairs. Notice that 
the results so far obtained are known and are also restricted to the standard 
statistics. 

Finally, the propagator <^ cLfc|,4T ^w'' 



I 

27r/ w 
i \ A 
2^ 



A 



2-7I-A2 
1 



2E, 



1 



w- Ek w + Ek 



(19) 
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where Ek — y{€k — IJ,)'^ + A| is the energy of the quasi-particle. When A = 

0, we get Ek = Sk — fJ', which is the electron energy with respect to the 
chemical potential fi and, indeed, superconducting solutions are given by A 7^ 
0. Furthermore, from Ek, one can easily realize that the minimum amount of 
energy to create an excitation in the superconducting state is A. 

Combining (16) and making use of the {w — w' i: iO"*")"^ — P {w — w')~^ ^ 
inS {w — w'), one have 



Z„J ^ '^'^ ' 2Ek J e(9-i)^^"' + l 



'^EkZc 



(20) 



One can easily realize that the whole g-dependence of the quantity above is 
contained in the generalized version of the Fermi distribution function, fq, 

since Z~^ Kq{^)Zi {^{q — 1)C/^) = 1 and fq can be also written in terms 
of its own integral representation. 

Finally, the self-consistent equation for the superconducting gap equation is: 

A 



^ = -Y.ykk'{bl)q=Y.Vkk' 

k k 



^Ek 



[1 - 2fq{(3Ek)] , 



(21) 



which is identical to that found by simply diagonalize the reduced Hamiltonian 
by a canonical variable transformation proposed by Bogoljubov[16] but with 
fq in the place of the usual Fermi function /. 

In the case usual superconductors it is reasonable to assume that the q involved 
in numerical calculations will not deviate appreciably from the unity, since 
most of their properties have fairly agreed with the BCS results in the standard 
statistics. Hence we may use in numerical calculations an approximated Fermi 
function[12]. 




where = [1 + (1 — q)x]~ is the generalized g-exponential function. The 
distribution function above avoids the self consistency found in (16) and might 
be considered a useful tool for systems near the extensivity[13]. 
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Following BCS, we replace the sum over k by an integral over 7 and the 
density of states shall be simply written as iV(7) ^ N{eF) or A^(0). Hence, 
the generalized self-consistency condition is 



hw]j 



N{eF)V 



and one can estimate the gap energy in the weak-coupling limit assuming 



ksTc 



2e-^« , (24) 



where 



~ d 



Ng^2 J dx -^fqi^) , (25) 



since [14]A(0) = 2hwD e-i/^(^F)v^ j^^ ^j^e limit g ^ 1 we get A(0)/A;bTc = 1.76, 
as usual. If one calculate self-consistently the gap equation taking into account 
the experimental data for Tc and 6^, one get that the coupling N{ep)V is much 
less than unity. Superconductors with larger ratio Tc/9d tend to have greater 
disagreement with the predictions of the BCS model in the usual statistics. 



4 Compcirison with Experimental Results 



In general, usual superconductors exhibit exponential behavior in the low tem- 
perature regime for several macroscopic measurable functions. As we will see 
below, the BCS model with its constant energy gap for a s-wave supercon- 
ductor (/ = 0) display such behavior in the standard statistics. On the other 
hand, the so-called exotic superconductors present power law behavior and it 
is believed that different pairing states may be responsible for such behavior. 
One of the most interesting results of the present paper is that, even for s-wave 
pairing state, it is possible to find power law behavior for different functions 
in the T <^ regime once q ^ 1- 

In this section we compare experimental results for specific heat, ultrasonic 
attenuation and tunneling experince for some usual superconductors with dif- 
ferent theoretical predictions for g = 1. We show that small deviations with 
respect to g = 1 provide better agreement when confronted to data, even 
regarding the approximations involved in the calculations. Furthermore, we 
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expect that the same q may fit the avaiable data for any different experiments 
of a given material. 



4-1 Specific Heat 



In an unitary volume, the eletronic specific heat in the superconductor state 
is 



(26) 



which may be interpreted as the specific heat associated to the quasi-particles 
with energy Ek, added to the contribution of the temperature dependent gap. 
Notice that when q 1, —-^fiix) — fi{x) {1 — fi{x)) and one recover the 
usual expression for a superconductor in the BG statistics. 

The "normalized" jump ( C| — C" )/C^ may be measured and compared with 
theoretical predictions. The BCS theory with standard statistics predicts the 
result of 1.43 for the specific heat jump, and usual superconductors exhibit 
some deviations from this value as it is shown on Table 4.1 for several ma- 
terials. We calculated numerically (C| — C*" )/Cg as a function of q and 
therefore we can estimate the pertinent value of q for the superconductors 
listed. Notice that was considered as a constant parameter in order to give 
the value of 1.43 at g = 1 for the specific heat jump. Notice, from Table 4.1, 
that the appropriate g-- value for Sn is 0.995. we shall see below that numerical 
calculations for Sn using q ~ 0.99 show better agreement with measurements 
for the specific heat in the low temperature regime, ultrasonic attenuation 
rate and tunneling experiences. The fact that the same q is able to give, at 
least qualitatively, a better agreement with different measurements serves as 
motivation to further application of nonextensive statistics in order to explain 
the anomalous features presented in some superconductors. 

Thus, according to the above discussion, we show the specific heat for Sn and 
V in Fig. 5. From Table 4.1, we expected that q = 0.995 should reproduce 
better data for Sn than q = 1. We see also that the specific heat measure- 
ments for vanadium and tin [26] are well fitted by the exponential function 
9.17e^^'^^^'"/^^ in the temperature region 0.25 < T/T^. < 0.7, which is quite re- 
markable if compared to the prediction of the usual BCS model, S.Se"^-^^^^^/^^ 
for 2.5 < Tc/T < 6. Theoretical plots ior q — 1 and 0.99 agree well to the 
experimental fitting above. However, one can see in Fig. 1, specially in the 
insert, that the agreement for V with q = 0.99 is better not only at T ~ Tc, 
but specially for lower temperatures. At such low temperatures there is an 
increasing departure from the standard statistics. 
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Table 1 

Display of the values q found when numerical plot is compared to the specific heat 
jump measured experimentally [22]. ci*'' was calculated numerically in the weak- 
coupling limit. Notice that these usual superconductors are well described by g 1. 

By the same token, on Fig. 2 we show several data for the specific heat for 
Al and Zn. They both deviate from the usual BCS calculations in the low- 
temperature regime. The data of Phillips[25] for Zn is better fitted by a T^- 
law in the low temperature regime. According to our estimations, shown in 
Table 4.1, the Al and Zn data should be fitted hj q = 1.03. Indeed, one can see 
in Table 4.1 that the calculations with q = 1.03 agree much better with the 
experimental results of several groups than the usual statistics. We can also 
see in Fig. 2 the tendency of the specific heat to deviate from a exponent to 
a powder law behavior even with such small value of g = 1.03. This behavior is 
even stronger with a much larger value of q, like q = 1.2 for instance, although 
such calculations will be presented elsewhere. 

4-2 Ultrasonic Attenuation Rate 

In this section we treat the problem of acoustic waves perturbing the electronic 
system described by the diagonalized BCS reduced Hamiltonian. We will con- 
sider only longitudinal waves and the temperature interval, the frequency and 
wave vector are related by the inequality w < vpp < A/h. As usual, we assume 
that the matrix element of the perturbation depends on k' — k. In a typical 
ultrasound experience the frequency is less than 10^ Hz, hence hw <tC ksT, 
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and in the limit of very low frequencies one get that the normalized ultrasonic 
attenuation rate within Tsallis' statistics is 



Notice that a. s join also has exponential decrease when — > 1 in the low tem- 
perature regime. Such behavior is experimentally verified in several usual su- 
perconductors, although plots of experimental data tend to fall off somewhat 
more rapidly just below Tg. Data of Morse and Bohm[29]for Sn are a example, 
Fig. 3, and are compared with theoretical plots for q — 0.99 and 1. Again, 
q = 0.99 is in better agreement with experimental results. An interpretation 
for the discrepancy between the usual theory and the data was given in terms 
of the gap function, which, compared with it value at T = 0, rises slightly too 
rapidly as the temperature decreases just below T^. Tunneling measurements 
for Sn by Morse and Bohm[30] are consistent with such feature. On Fig. 4 we 
show their results for Sn and In with the calculation with q — 0.99 and 1. 
We can easily see that q = 0.99 is in much better agreement than the usual 
q = I results. Thus in the present interpretation this result is attributed to 
the nonextensivity of the system. 



5 Conclusions 

We conclude this paper pointing out our new findings: 

• We have show that the BCS potential widely used in the description of 
superconductors is nonextensive. Following this result we have developed 
a BCS theory with nonextensive entropy which depends on the index q. 
We have, therefore, derived some thermodynamic quantities that can be 
compared with the experiments. 

• Using the experimental values of the specific heat, we estimated the nonex- 
tensive index q for several superconductors. Explicit calculations with such 
values of q for several materials arc indeed in much better agreement with 
the present experimental results than the usual values. In particular for Sn 
we have estimated q ~ 0.99 and have shown that with such value we re- 
produced with a better agreement the experimental results of the specific 
heat, ultrasonic attenuation and tunneling experiments. Similar calculations 
were performed for Al and Zn. The fact that three well-known properties 
are fitted by the same q 7^ I rules out a possible deviation due to the 
approximations on the calculations. 




(27) 



—00 
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Furthermore, we have seen that as q increases, the exponential behavior of the 
specific heat changes in a power law behavior. This is a novel result and it 
means that even if there is a constant gap in the neighborhood of a Fermi sur- 
face, we can get a power law behavior once we use Tsallis statistics with q ^ 1. 
Several new superconductors, like high temperature superconductors for in- 
stance, exhibit such power law behavior which is attributed to a non-constant 
d-wave order parameter. At the moment we are developing a non weak-coupled 
theory of superconductivity with Tsallis statistics to deal with these and oth- 
ers exotic superconductors with thermodynamical quantities which exhibit a 
power law behavior. 
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Fig. 1. Comparison between the fitting of the low temperature electronic specific 

heat of Tin, ci"^ = 9.17exp[ -1.5( T/Tc ) ] ['yTc], and the theoretical plots for g = 1 
and 0.99. Both values of q describe very well the experimental results, but for V 
in particular, regarding the interval of temperature 0.25 < T/Tc < 0.7, a better 
agreement was found with q = 0.99. 



Fig. 2. Experimental data for the specific heat of Al and Zn [23]. Notice that 
the measurements by Phillips[25], Goodman[24] e Zavaritskii[27] for the Aluminum 
are well described by the generalized BCS model with q = 1.03, the same value 
furnished in Table 4.1. Also notice that data obtained by Phillips[25] for Zinc are 
very well described by 3{Tc/T)^ in the low temperature regime. 



Fig. 3. Ultrasonic attenuation rate calculated for different values of q and compared 
to experimental data for Sn[29]. A better agreement was found with q = 0.99, as 
proposed by measurements of the specific heat jump. 



Fig. 4. Comparison between the gap energy from tunneling measurements for Tin 
and Indium[30] and theoretical plots using q = 1 and 0.99. Notice that the gap 
function related to q = 0.99 have better qualitative agreement with experimental 
results. 
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